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Abstract. We develop a model of thin turbulent accretion discs supported by magnetic pressure of turbulent 
magnetic fields. This applies when the turbulent kinetic and magnetic energy densities are greater than the 
thermal energy density in the disc. Whether such discs survive in nature or not remains to be determined, but 
here we simply demonstrate that self-consistent solutions exist when the a-prescription for the viscous stress, 
similar to that of the original Shakura-Sunyaev model, is used. We show that a ~ 1 for the strongly magnetized 
case and we calculate the radial structure and emission spectra from the disc in the regime when it is optically 
thick. Strongly magnetized optically thick discs can apply to the full range of disc radii for objects 10 -2 of the 
Eddington luminosity or for the outer parts of discs in higher luminosity sources. In the limit that the magnetic 
pressure is equal to the thermal or radiation pressure, our strongly magnetized disc model transforms into the 
Shakura-Sunyaev model with a = 1. Our model produces spectra quite similar to those of standard Shakura- 
Sunyaev models. In our comparative study, we also discovered a small discrepancy in the spectral calculations of 
Shakura & Sunyaev (1973). 
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1. Introduction 



The well known and most widely used model of the accretion disc was proposed and elaborated b vlShakural h 972ft and 
Shakura fc Sunvaevl l) 19 73ft . In this model the disc is vertically supported by the thermal pressure l|Shakura fc Sunvaevl 



1973ft . Turbulent viscosity is invoked in the Shakura-Sunya ev model to explai n the angular momentum trans fer 



required by the accretion flow. As originally pointed out in iLvnden-B cll (1969J) and IShakura fc Sunvaevl l|1973ft a 
magnetic field can also contribute to the angular momentum transport. A robust mechanism of the excitation of 
magnetohydrodyn amical (MHD) turbulen ce was shown to operate in accretion discs due to the magneto-rotational 
(MRI) instability l)Balbus fc Haw lcv. 1998). The growth of the MRI leads to the excitation of turbulent magnetic fields 
and self-sustained MHD turbulence. The contribution of Maxwell stresses to the transport of angular momentum is 
usually larger than Reynolds stresses. However, the ma gnetic energy observ ed in many numerical experiments was 
smaller than the thermal energy of the gas in the disc ijBrandenburd. Il998ft . Simulations of the non-linear stage of 
MRI are typically local simulations in a shearing box of an initially uniform small part of the disc. Attempts to 
expand the computational domain to include a w ider area of radii and azimuthal angle llHawlev fc KrolM l200lt 
Hawlev, 2001; Armitagc, Reynolds & Chiang, 2001) are underway. However, even before the recent focus on the MRI 



Shibata. Taiima fc Matsumotol jl990ft observed the formation of transient low j3 state in a shearing box simulations 



of the non-linear Parker instability in an accretion disc. 

Vertical stratifica tion has been considered in the s hearing box approximation l|Brandenburg et "all [l995; 
iMiller fc Stonel 12000ft . In particular. iMiller fc Stond l)2000ft investigated discs with initial Gaussian density profiles 
supported by thermal pressure. The initial seed magnetic field grows and starts to contribu te to the vertica l press ure 
balance. The computational domain extends over enough vertical scale heights to enable IMiller fc Stonel l|2000ft to 
simulate the development of a magnetically dominated corona above the disc surface. In the case of an initial axial 
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magnetic field, iMiller fc Ston3 (|2000^ observed that the saturated magnetic pressure dominates thermal pressure not 
only in the corona but everywhere in the disc. As a consequence, the thickness of the disc increases until it reaches 
the axial boundaries of the co mputational box. The formation of low 8 filaments in magnetized tori was also observed 
in global MHD simulations by iMachida. Havashi &: Matsumotol l)2000|) . Although further global MHD simulations of 
vertically stratified accretion discs are needed, this numerical evidence suggests that magnetically dominated thin discs 
may exist. 

Previously , analytic models of thin acc retio n discs with angular mome ntum transfer due to magnetic stresses were 
considered bv lEardlev &: Lightmanl l)l975(l and lField k, Rogers! lll993afb!h Both these works included magnetic loops 
with size of the order of the disc thickness. In lEardlev fc L ightmanT I|l975(K the magnetic loops were confined to the 
disc. Loop stretching by differential rotation was balanced by reconnection. The reconnection speed was a fraction 
of the Alfven speed. Radial m agnetic flux was cons i dered as a free function of the radius. Vertical equilibrium and 
heat transfer were treated as in lShakura k Sunvaevl (^973), with the addition of the magnetic pressure in the ver tical 
support. No self-consistent magnetically dominated solutions were found in model of lEardlev k Lightmanl ljl975|) . 

In contrast, d ominance of the magneti c pressure over the thermal and radiation pressure was postulated from 
the beginning by iField k Rogers! l(l993alrj l and verified at the end of their work. These authors assumed that the 
ordered magnetic field in the disc, amplified by differential rotation, emerges as loops above the surface of the disc 
due to Parker instability. Because the radial magnetic field in the disc has an intially sectorial structure, the loops 
above th e disc come to close cont act and reconnect. All dissipation of magnetic field occurs in the corona in the 
model of IField fc Rogers! lll993albh . Such a corona was assumed to be consisting of electrons and some fraction of 
positrons and no outflow from the disc is present. Electrons and positrons are accelerated to relativistic energies at 
the reconnection sites in the disc corona and su bsequently emit synchrotr on and inverse Compton photons. Because 
reconnection was assumed to occur at loop tops, IField fc R.ogersl (|l993alb| ) found that up to 70 per cent of the energy 
released in reconnection events in the corona will be deposited back to the surface of the disc in the form of relativistic 
particles and radiation. Only thin surface of optically thick disc is heated and cools by the thermal emission, which is 
the primary source of soft photons for the inverse Compton scattering by relativistic particles in the corona. 

Since the characteristic velocity of rise of the loops of the buoyant magnetic field is of the order of the Alfven speed , 
it takes about the time of a Keplerian revolution for the loop of the magnetic field to rise fe.g.. lBeloborodovl [l999'). 
This is also about the characteristic dissipation time of the mag n etic field in shocks inside the disc (see section |2J). 
The model we explore here differs from that of IField k Rogers! l|l993alb| ) in that the dissipation of the magnetic 
energy occurs essentially inside the disc and the heat produced is transported to the disc surface and radiated away. 
Observations of hard X-ray flux indicate the presence of hot coronae where a significant fraction of the total accretion 
power is dissipated. For example, the X-ray band carries a significant fraction of the total luminosity of Seyfert nuclei: 
the flux in the 1—10 keV band is about 1 /6 of the total flux from IR to X-rays, and the flux in 1-500 keV band is about 
30-40 per cent of the total energy output ijMushotzkv. Done fc Pounds! 1 1993|K Another example is the low/hard state 
of galactic black hole sources, where the borad band spectrum is completely dominated by a hard X-ray power law, 
rolling over at energies of ~ 150 keV l|Nowakl |l995t IrJoneL 120021 . Also, in the so called very high state, some of galactic 
black hole X-ray sources show both t hermal and no n-thermal (power law) components, with the ratio of non-thermal 
to total luminosity of 20-40 per cent l|Nowaklll99l . Reconnection events and particle acceleration should also happen 
in rarefied strongly magnetized corona of the disc in our model and could cause observed X-ray flaring events. However, 
we do not consider the coronal dynamics here, and instead just focus on the structure and the emission spectrum of 
the disc itself. 

Models of magnetized accretion discs with externally imposed large scale vertical magnetic field and anomalous 
magnetic field d i ffusion due to enhance d turbulent diffusion have also been considered ( Shalvbkov fc Rtidigeri l2000l: 
ICamobellL EooH iQgilvie fc Livid. l200lh . The magnetic field in these models was strong enough to be dynamically 
important. But those models are limited to the subsonic turbulence in the disc and the viscosity and magnetic 
diffusivity are due to hydrodynamic turbulence. Angular momentum transport in those models are due to the large 
scale global magnetic fields. Both small scale and large scale magnetic fields should be present in real accretion discs. 
Here we consider the possibility that the magnetic field has dominant small scale component, that is magnetic field 
inside the disc consists mostly of loops with size less than or comparable to the thickness of the disc. 

We consider vertically integrated equations describing the radial structure of the magnetically dominated turbulent 
accretion disc and provide the solutions for the radial dependences of the averaged quantities in section[2| In section|2] 
we analyse the conditions for a magnetically supported disc to be self-consistent. In section 01 we calculate thermal 
emission spectra of magnetically supported disc taking into account scattering by free electrons. 



2. Radial Disc Structure 



Here we neglect effects of general relativity and do not consider the behaviour of the material closer than the radius 
of the innermost circular stable orbit r s . We assume a non- rotating black hole with r s — 3r g , where the gravitational 
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radius of the black hole of mass M = Mg x 10 8 Mq is r g = 2GM/c 2 — 3 x 10 13 cmM 8 . We assume that accretion 
occurs in the form of a geometrically thin accretion disc and verify this assumption in section [21 We consider a disc of 
half-thickness H, surface density S, averaged over the disc thickness volume density p = £/(2iJ), accretion rate M, 
and radial inflow velocity v r , v r > for the accretion. We take fl K = (GM/r 3 ) 1 ' 2 = 7.2 x 1CT 4 s~ 1 Mg 1 (r g /r) 3 / 2 to 
be the angular Keplerian frequency at the radial distance r from the black hole. Then, equation of mass conservation 
reads 

M = 2wrT,Vr. (1) 

In the stationary state M does not depend on time. Equation of angular momentum conservation is (e.g., 
IShapiro fc Teukolskvl Jl983^ 

f^{r)2H2TTr 2 = M (VGMr - (y/GM^j , (2) 

where is the tangential stress at a radius r, which acts on the inner part of the disc, factor £ accounts for the 
unknown torque acting on the disc at the inner edge r s . In the standard Shakura-Sunyaev model, £ — 1, which 
corresponds to a zero torque inner boundary. The tangential stress at the innermost stable orbit, /^(?" s ) is related to 
the parameter £ as follows 



U(r s )2H2Trr 2 - (1 - QMy/GMr,. (3) 

The stress fcf,(r s ) > 0, if the torque spins up the gas at r s , and f<j>(r s ) < 0, if the torque retards the rotation of the gas. 
Since the magnetic and turbulent pressures exceed the thermal pressure, by assumption, we neglect thermal pressure 
contribution to the vertical pressure balance. Assuming equipartition between turbulent and ma gnetic pressure in 
supersonic MHD turbulence l|Stone. Ostriker fc Ganinfiel Il998t IStoneL Il999t iMiller & Stcmel 1200(11) it is easy to show 
that the equation of vertical equilibrium is solved to give approximately the result 



H K Ok V 2ttS ' (4) 

where va — -, tt-t^ is the average Alfven velocity in the disc. We use the 'a-prescription' for the magnetic viscosity 

(47rp) 1 ' 2 

B 2 

in the disc, i.e. taking the tangential stress to be proportional to the sum of the magnetic pressure, — , and the 

8ir 

turbulent pressure, pv 2 /2. We do not consider global transport of angular momentum and energy across the disc due 
from large scale magnetic fields. We assume the magnetic field to be sufficiently tangled on scales of order of the 
thickness of the disc, such that large scale field is small compared to this tangled field produced by turbulence. In 
B 2 

equipartition, — w pvz/2 and a-prescription becomes 

87T 

U = af^- (5) 

The dissipation of the magnetic and kinetic energy causes heat input in the disc which is balanced by the heat losses 
due to radiation. If the cooling is efficient enough such that the time of radial advection of the heat due to the accretion 
flow is much longer than a few Keplerian periods, the heating and cooling balance, and establish an equilibrium kinetic 
temperature in the disc. However, we assume that this temperature is insufficient for the associated thermal pressure 
to contribute significant vertical support. The total radiated energy from the unit surface of the disc will be the same 
as in the standard Shakura-Sunyaev model. This energy flux i s independent o f the viscosity mechanism assumed, but 
depends on the inner torque boundary condition (see also iGammiel ()2000|) for the effects of varying the torque at 
r s ). Thus, the system of equations JIJ, J2|, (@J and JSJ decouples from the energy balance equation. 

Equation (0J can be solved to give H as 



H =2^- ^ 
With the prescription for the viscous stress, equation JSJ), the angular momentum conservation equation J5J) becomes 



It is remarkable that the value of a can be constrained in our model. In the framework of the local approach all 
the work done by non-gravitational forces on a patch of the disc is reduced to the work done by viscous stress f^. 
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Physically, this work results from the action of turbulent and Maxwell stresses. Heating occurs from dissipation of 
supersonic turbulence. The rate of such heating can be expressed through the kinematic viscosity coefficient v in the 
usual way (heating per unit volume) 

PT% = f X (8) 
at vp 

where T is the average temperature of gas inside the disc, s is the entropy per unit mass (e.g., see in 
IShapiro fe Teukolskvl l|l983|) l. Viscous stress, for Keplerian shear is 

3 

U = (9) 
Comparing this to the alternative expression (J5J for we see that a and v are related by 

v =«ir/- (10) 

Using expressions © and IjlOjl to substitute in equation JBJ) one can obtain the rate of heating per unit area of the 
disc Q as follows 

ds 

Q = 2H P T -r = 3apv 3 A . (11) 
One can also rewrite Q as 

y 

B 1 o o 

density of the turbulent pulsations and magnetic fields is — h —pv t — pv A under the assumption of cquipartition. 

In the stationary state, this turbulent energy dissipates with the rate Q. Therefore, the characteristic time of the 
dissipation of the turbulence is 

- -Q- ~ 3?w (13) 

On the other hand, dissipation of the supersonic turbulence occurs on the time-scale r t of the flow crossing the 
largest flow coherence size It of the turbulence. Direct dissipation at the shock fronts dominates the turbulent cascade 
of energy down to the microscopic resistive scale and leads to the enhanced rate of the dissipation. The question of the 
dissipation of supersonic MHD turbulence has been studied in connection with the interstellar turbulence, which is 
observed to be highly supersonic. Direct numerical simulations of both stea dy-state driven and freely decaying MHD 
turbulence l|Stone et alJll998t Stonelll999l lOstriker. Stone fe Gamniie]l200l|) all confirm this picture. Even if initially 



4ir r 3 \ V r 



the same expression familiar in the standard disc model (e.g., see IShapiro fc Teukolskvl £983)). The total energ; 



the motion is set up to be incompressible in the numerical experiments, shocks develop rapidly due to the non-linear 
conversion of the waves and the dissipation becomes dominated by the dissipation on shocks. This dissipation time 
is r t = h/vt- The coefficient of turbulent viscosity v is also related to the largest flow coherence size and turbulent 
velocity as 

v = ^kvt- 

We already know that v t — i>a from equipartition of magnetic and kinetic energies in turbulence. Also, v is expressed 
through a by relation (|10|l . This allows to estimate the largest flow coherence size of the turbulence 

l t = 2aH. (14) 

2a 

Consequently, the flow crossing time of the largest coherence size becomes Tt = . It should be that r t = ta with 

"K 

ta given by expression H13[). This is possible only when a = 1/V3- Due to approximate nature of all calculations 
which lead us to this value for a, it is only meaningful to state that a should be of order of 1. We assume a = 1 for 
all further estimates. For a w 1 the largest flow coherence size of the turbulence becomes of order of the thickness 
of the disc It ~ 2H, and the turbulent viscosity coefficient takes its largest possible value v ps Hva, still compatible 
with the local viscous description of the disc. The dissipation time-scale for the magnetic turbulence is ta ~ It 
is very probable that such large scale turbulence will lead to the buoyant rising of the magnetic field loops into the 
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corona, subsequent shearing by the differential rotation and reconnecti on of the loops. This can r esult in the formation 
of the hot corona or acceleration of particles to relativistic energies (iField fc Rogers! Il993albfr . The formation of a 
magnetized corona and the emission spectrum from the corona are important, however, here we focus on the disc. 

The free parameters are M and £. Also, we need to specify one more physical condition, since the dependence of B 
on radius in equation Q is undetermined. Such a condition should come from the modelling of supersonic turbulence. 
Lacking a detailed model, we assume that the radial dependence of the vertically averaged magnetic field in the disc 
is the power law 



-<5 



B = Bw {^k) ' (15) 

where Bio is the strength of the magnetic field at 10r g , S > is some constant. Accretion rate M can be related to 
the total luminosity of the disc L and the radiated fraction e of the rest mass accretion flux Mc 2 . At the luminosity 
of an AGN 

L = 1.3 x 10 46 Z E M 8 erg s"\ (16) 
the mass flux is 

M = 0.23 MQyr" 1 f M 8 = 

l^xlO^gs-^^Mg. (17) 

Here we denote the ratio of disc luminosity to the Eddington luminosity by = L/L c ad- The value of e is determined 
by the inner boundary condition at r = r s . Typically, e « 0.1. Using such parametrization and the expression 115fl for 
the magnetic field, one can obtain the following radial profiles of H, £ and p from equations 10 and Q 

h n K Aig x z E / b 10 y 2 

— = 2.1 x 10 — — 7—- x 



r g aB 2 r g ' 2e V 10 4 G 



/ \ 2(5-3/2 



V 10ro 



-l 



„ Q-B , o / Ir\ 

£ = — = 5.1 x 10 3 g cm" 2 — x 



B 



10 



9/2-4(5 



u.'g; M * Q ~\mJ ■ (19) 

a " B6 ■ -n-10„_-3 fl^" 2 



p = , = 4 x 10" 1U g cnr~ J — x 

/ l) \ 6 / \ 6—6(5 

' Bw ^ G~ 2 M% ( ^) , (20) 



V 10 4 G / V 10r, 



where Q = 1-CW^. We see that the disc becomes geometrically thicker and less dense when magnetic field decreases: 

weaker magnetic field leads to weaker tangential stress (equation JSJ); H increases in order to accomodate constant 
angular momentum flux for the same M such that H oc B^ 2 (equation (J5J); T, and p should decrease strongly, S oc Bf 
and p cx Bf , in order to ensure vertical equilibrium with larger H and less pressure support from B 2 (equation I©); 
radial inflow velocity increases as v r oc B^ 4 to ensure constant mass flux. 

Let us now summarize the similarities and differences between our model and the standard Shakura-Sunyaev model. 
If we replace the thermal pressure in the standard model by the sum of the magnetic and turbulent pressures, the 
equations for mass conservation l(T|l. angular momentum conservation l|2|l. the viscosity prescription and vertical 
pressure support 10} are the same as in the standard model. The pressure in the standard model is determined by the 
rate of the cooling of the disc, while the a coefficient can be an arbitrary function of r, a(r) < 1. In our model we have 
the magnetic pressure unspecified in its radial dependence as soon as it exceeds the thermal pressure, but a ~ 1 for 
all r. The latter results from the much faster dissipation of supersonic turbulence than subsonic turbulence assumed 
in the standard model. Both our model and the standard model have only one undetermined function of radius, (a(r) 
in the standard model and B(r) in our model). The determination of this free function would eventually come from 
detailed modelling of the MHD turbulence. 
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3. Estimates of the Disc Parameters 

Now let us use the solution for the disc structure provided by equations 0, 10 and (|15fl and obtain constraints on free 
parameters of the model, such that our model of thin magnetized accretion disc is self-consistent. Using equations 
and to substitute for M in equation (0 we can express the radial inflow velocity as 

- = ( 21 ) 

vk y r z 

The factor Q vanishes at r = r s for £ = 1 and so the surface densi t y E o f the disc diverges near r = r s . The same 
divergence also occurs in standard disc model (|Shakura & Sunvaevl 197.^ . In reality, of course, viscous torque does 
not vanish at r = r s , £ ^ 1 but close to 1 and there is no divergence of S at r = r s . Determining £ would require 
a full general relativistic treatment of the accretion flow close to the black hole and is beyond the immediate scope 
of this work. Only parts of the disc close to r s are sensitive to the exact value of C- For r > r s , the disc structure is 
approximated well with the formulae in section [2] for f = 1. Provided that r > r s and a < 1 one can see that the 
radial inflow velocity for a thin disc (H <C r) is always a small fraction of the Alfven velocity which, in turn, is a 
small fraction of the Keplerian velocity. Therefore, the dissipation time-scale ta is always much shorter than the radial 
inflow time-scale. In the stationary case this means that the advective terms in the energy balance equation can be 
always neglected. Energy balance becomes local: the rate of gas heating Q should be balanced by the cooling due to 
radiative processes. Now we consider the physics of radiative cooling which determines the disc temperature. 

A necessary condition for the existence of a magnetically dominated disc is that the vertical escape time for 
radiation must be shorter than ta, so that the energy density of radiation, e r = aT 4 , remains smaller than the energy 

B 2 

density of the magnetic field £b — — ■ For optically thin, geometrically thin discs this condition is always satisfied 

87T 

since the inverse of escape time c/ H 3> ^k- As we will see, Thomson scattering is the dominant source of opacity in 
most cases of optically thick discs. The average time it takes for a photon to escape out of optically thick disc with 
optical depth r 3> 1 is tH/c. For Thomson scattering r = t c = Hn e aT = ktE/2, where n e is the number density of 
free electrons, or is the Thomson cross section and kt = 0.4cm 2 g _1 is the Thomson scattering opacity. For simplicity 
we assume the composition of the disc to be completely ionized hydrogen. Then, n c is equal to the number density of 
protons in the disc, n. The necessary condition now becomes 

H 2 1 

na T < — . (22) 

C ilK 

Using equations (0 and one can rewrite the condition ll'L'l) as 

< I- (23) 

We express this and all the subsequent conditions in terms of free parameters of the model: Bio, 8, 1-e/c, and Mg. 
Using the expression l|18|l for H, the necessary condition IL'-it becomes 

/ D \ 2 / \ 3/2 — 28 

"* {wk) "> e iik) <<L (24) 

The condition for the disc to be optically thick for Thomson scattering r c ^> 1 becomes 
2.0xl0 3 (^y 1 f^Vx 



2e) V!0 4 G 

/ \ 9/2-45 

M ^ _1 [ik) ^ L (25) 

In an optically thick disc, radiation is transported by turbulent motions and radiative diffusion. The characteristic 
time-scale for the turbulent transport of radiation to the surface of the disc is ~ 1/^k- When condition (|22() is 
satisfied, the diffusion of radiation dominates the advection due to turbulent motions. Thus, we can neglect turbulent 
convective transport of radiation for any values of parameters, whenever steady state magnetically dominated model 
of the accretion disc is considered. In the steady state, the radiation flux from each surface of the disc must be equal to 
Q/2 with the dissipation rate Q given by expressions or Q12[l. The effective temperature of the escaping radiation 
flux is determined by this dissipation rate Q as 

T„ - (f) V4 - 7.5 x to* (£)■",< 

(is;) ff,,4M » 1/4 - (26 » 
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and is the same in our model as in the standard disc model. True absorption of photons in free-free transitions also 
occur in the disc. For an approximate estimate of the radiative conditions in the disc we consider Rosseland mean of 
the free-free absorption opacity 

% = 6.4 x 10 22 cm 2 g- 1 p7 1 - 7/2 , (27) 

_o O 

where p is expressed in g cm 6 and T is expressed in K. When the effective optical thickness of the disc 2t* = 
S-^/ (act + ks)ks is large, then local thermal equilibrium is established in the disc and the radiative flux is described 
by diffusion approximation. Only in the thin surface layer at a distance from the surface less then the thermalization 
optical depth f* » 1, the spectrum of radiation deviates substantially from that of the black body (see section0]below) . 
The solution of the diffusive radiation transport gives the usual result relating the temperature at the midplane of the 
disc, T mp( j, with the effective temperature at the surface llShakura fc Sunvaevt Il973t iKrolikl fl999^ 

T mpd « T cS (^^] 1 = 3.3 x 10 5 K x 



16 / 



Bin \ ,.1/4/ r 



3/8-5 



io<g, m ' UsJ ■ < 28 > 

The averaged temperature in the disc is close to T mp d with the actual profile being determined by the details of the 
vertical dependence of the dissipation rate in shocks. As in the case of the inner radiation dominated part of the 
standard disc, T mpc j does not depend on the accretion rate. However, it is directly proportional to the value of the 
magnetic field and has radial dependence governed by the 5. 

The dominance of the magnetic and turbulent energy compared to the energy density of radiation is expressed as 
a ^m P d ^ P V A- O ne can substitute here for T mp d from equations (|2*%)l and J2HJ. Heating rate Q is given by equation (fTTJl . 
After using equations H18H20[1 to manipulate with S, B and H, one can reduce the condition of magnetic pressure 
dominance over the radiation pressure to be 

9 * nH i«-L. (29) 
4c SZk 

This differs only by a factor of order unity from the necessary condition for the escape of radiation, equation H22[). The 
assertion that the condition (|24|l implies smallness of radiation pressure is thus confirmed by direct calculation. The 
next condition we need to impose is that the gas pressure is small compared to the magnetic pressure and turbulent 
stresses, which is equivalent to the statement that the turbulence is highly supersonic. For a thermalized radiation 
field, this condition is nfcT mp d <C pv\. Using expressions (|28|l and (|26|l for the temperature and substituting for the 
density and H from expressions <|18H2Uf) one can express the condition nfcT mp d >C pv\ as 



4_ \ V2 

a 9 / 2 n^ 17/2 B 10 M- 4 g- 4 « 1. (30) 



32ir 2 acm 5 



p ■ 

Substituting for fix, B, and M in 13U|) and taking the —1/5 power we obtain 



-51/20+25 



g 



Because the free- free opacity Rff strongly depends on temperature an d density (equat ion (l2*7|l). the vertical density 
and temperature distributions are needed to evaluate f*. Simulations in iMiller fc Stone] ll200fll) show a s harp density 



drop off by two orders of magnitude at the surfaces of the disc (see figure 11 in IMiller fc Stonel 1)200(1) ) . Inside the 
slab bounded by this density drop off the magnetic field and kinetic energy are in approximate equipartition. For the 
purposes of calculating effective optical depth for absorption we approximate the density profile between z = —H 
and z = H inside the disc as constant and assume zero density at the disc surfaces at z = ±H . Modest variations of 
magnetic field and density across the disc height observed in numerical simulations support this and also suggest the 
assumption of a uniform turbulence dissipation rate across the thickne ss. With these approximations, the soluti on for 
diffusive radiation transport in the vertical direction z can be written ijShakura fc Sunvaevlll973llKrolild . 17999) 



rp4 rp4 

mpd 



Z 

3r c ' V~H 



(32) 
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The Eddington approximation near the surface z = H of the disc was used to obtain this temperature distribution, 
and T e ff = \[2 T\ Z _ H according to the solution l|32l) . This solution is not valid near the surface of the disc in the region 
dominated by Compton scattering but gives T in the bulk of the disc if f* ^ 1. Using the expression l|32[l for T one 
has for the optical depth of the disc for free-free absorption 



r ff = / P Kg dz = 6.4 x 10 22 P 2 T- 7/ d 2 H 
Jo 



1 / on -7/8 



i - £ + ) ( 33 ) 

where £ = z/H, p and T are expressed in g cm -3 and K. The value of the integral over £ determines how much the 
actual value fg is larger than the value of fg obtained if one assumes T — T mpc i everywhere in the disc. The £ integral 
in equation is calculated to be 

1 / o \ -7/8 

_8_y /8 _(_?_ xl/8 

3t c / V 3r ( 

Because r c 3> 1 the number in square parenthesis is close to 1, so we can omit it and obtain the final result for fg 

f ff = 5 x 10 2 V^f (34) 

When one substitutes for T mp( j, H, and p their expressions through the magnetic field and the accretion rate, the 
expression (|34[) becomes 

147/16-135/2 

f ff = 8.6x 1(T 2 M 25/8 f ^) G- 3 



10r e 



j \ -3 / R X 13/2 



2e7 \10 4 G 



(35) 



Generally for 5 w 1, Mg ^ 1, ZeA ~ 0.1, -Bio 10 4 G, and r ~ 10r g , fg ~ 1. However, because of steep dependence 
of fg on Ie, Ms, and Bi , the value of Tg can become large for lower accretion rates, more massive black holes, and 
stronger magnetic fields. The ratio of fg to r c is 

/ x 75/16-5(5/2 

^ = 4.2xl0- 5 M 8 9/8 (^-) <T 2 
Tc. V lOrv 



'E \ / i>10 



2eJ V10 4 G 



(36) 



We see that for the typical values of the parameters above the ratio fg/r c ~ 10 -2 but free- free optical depth can 
exceed Compton scattering optical depth for smaller values of Ie and at larger radii. When fg <C r c the effective 
optical thickness of the disc becomes f* = \/fgr c or 

219/32-21(5/4 

is x m 8 41/16 ( — ) g x 

Ze\ / -Bio 



2eJ \10 4 G 



10r g 

21/4 

(37) 



When Tg ^> t c the effective optical thickness is equal to Tg. 

Finally, the ratio of the disc semi-thickness H to the radius r using (|18fl is 

^ = 2.1 x 10- 2 ^ ( JQ) 2 M-'G ( 2 ^ 5/2 . (38) 
r 2e Vl0 4 Gy 8 \lOr s ) K 1 

Now we summarize conditions when our model is valid: 



1. The necessary condition (|24|) . which is also the condition for the dominance of the magnetic pressure over the 
radiation pressure. 
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2. The condition (|31|) for the dominance of the magnetic pressure over the thermal pressure of gas in thermalized 
optically thick disc. 

3. The condition for the disc to be optically thick. This is either condition l|25|l t c ^> 1 or the condition that fg given 
by expression l|35[> is greater than 1. 

4. The condition that the gas and radiation inside the disc are in local thermal equilibrium, r* 3> 1, where f„ = fg if 
Vs l T c > lj an d f* is given by expression (|37[1 if fg /r c < 1 . 

5. H/r < 1. 

We varied the parameters Ze/c, S, Ms to obtain the allowed region of our disc model in the -Bio, r / r g plane, using the 
above five conditions. These plots are shown in Figs. El 

Depending on the power S in the dependence of magnetic field B cx r~ A , optically thick magnetically dominated 
accretion discs can exist only at a limited interval of radii. For 6 ~ 1 (as on our plots for i5 = 5/4) a thin magnetically 
dominated disc (shadowed regions on the plots) is possible for 5r g < r < 1000r g . The window for the strength of 
the magnetic field is not very wide: about one order of magnitude or less. This window is narrower for low masses of 
the central black hole and is wider for higher masses. The value ZE/(2ae) = 10 corresponds to about the Eddington 
accretion rate for e 0.1, and because of a w 1 in our model (see section 2). Higher values of l-E/(2ae) correspond 
to higher accretion rates. Allowed values of the magnetic field are higher for higher accretion rates. The magnetic 
fields in the discs around higher mass black holes are smaller than in the discs around lower masses black holes as 
is temperature of the disc (T mp( j) and the surface radiation flux. For large luminosities (Ze/(2cv£) <; 1) the inner disc 
cannot be optically thick for true absorption but can be optically thick to free electron scattering. Comptonisation 
becomes significant for in the inner regions at high luminosities (see section^J. We leave consideration of Comptonised 
regimes for future work. 

Physically, the limitations on the magnetic field strength can be understood as follows: suppose one decreases 
Biq while keeping M constant. Then, H cx B^ is increasing; £ cx Bf , and p cx Bf , both decreasing (equations 
Scattering opacity through the disc r c cx T, strongly decreases, so the heat is transported to the surface faster 
and T mp( j cx Bio decreases (equation l)28p): thermal and radiation pressures decrease as Pth oc B\ and P ra( j cx Bf 
respectively; plasma parameter j3 decreases, so the disc becomes more magnetically dominated; fg and f* both decrease 
since their decrease due to lower E overcomes the increase of the absorptive opacity from the drop of the temperature. 
Therefore, there exists only a limited interval of Bio such that (3 < 1 still the disc is optically thick to true absorption. 



4. Radiation Spectra of Optically Thick Magnetically Dominated Disc 

4.1. Calculation of spectra 

Free-free, bound-free and cyclotron emission could contribute to the radiation spectrum. In Appendix A we show 
that, because of the self-absorption in the dense disc, the total flux of cyclotron emission from the disc surface is 
negligibly small compared to the total radiated power Q. This power is entirely due to free-free and bound-free 
radiative tr ansitions. Cyclotron and synchrotron emission can be importan t in the rarefied and strongly magnetized 
disc corona jlkhsanovUlQiS iField fc Rogerslll993aL iDi Matteo eraflll9974 , but our focus here is on the disc. 

We perform a simplified calculation of the emergent spectrum. We assume local thermodynamic equilibrium and 
do not consider effects of the temperature change with depth. This is justified when the spectrum is formed in thin 
layer near the disc surface. For simplicity we do not include the bound-free contribution to the opacity. Free-free 
opacities are the dominant source of thermal absorption for T 10 5 K, so our simplified spectrum is most relevant 
for smaller masses of the central black hole, for which the inner disc is hotter. Our goal here is to capture the effect 
of the magnetic field on the shape of the spectrum. We consider only optically thick disc models with both t c ^> 1 
and r* 3> 1. The electron scattering opacity kt does not depend on frequency in the non-relativistic limit, whereas 
free- free absorption opacity Kg(v) is a function of frequency: 

Ks (v) = 1.5 x 10 25 cm 2 g~Vr~ 7/2 1 ~ 3 g(x, T), (39) 

where we denote x — — . Thus, the parameters of the disc model Bio and 8 will affect the spectrum by means of Kg 
kT 

dependence on p and T/T e g. Gaunt factor g(x,T) is slowly varying function of a; and T between approximately 0.3 
and 5 in wide range of frequencies and temperatures l)Rvbicki fc Lightmanlll979l) . Moreover, in the wide temperature 
interval 10 3 K < T < 10 s K Gaunt factor is between 0.5 and 2 for frequencies 0.1 < x < 10 near the maximum 
of thermal emission. It is quite reasonable to set g(x,T) = 1, which we do in all further calculations. Note, that 
Kg behaves as 1/a; 2 for x <C 1, so free-free absorption is always more important at lower frequencies than electron 
scattering. 
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Fig. 1. Plots of conditions when our model is valid. Plots are for M = 10 8 Mq and 6 — 5/4. On each plot values of 
radius extend from 5r g to 1000r g . The panels differ by values of l-E/{2ae). On each panel, any given model of the disc 
is represented with a horizontal line B w ~ constant. Filled areas indicate regions where a magnetically dominated 
geometrically thin and optically thick disc can exist. The difference in filling represents different types of spectra 
emitted locally from the surface of the disc: the regions with black body spectra are filled with lines in top-left to 
bottom-right direction and the regions with modified black body spectra are filled with lines in bottom-left to top-right 
direction. There are seven types of lines on each panel: (1) upper solid line curved upward on each plot bounds the 
region where radiation pressure is small compared to magnetic pressure and turbulent stress (below the curve); (2) 
lower solid line curved downward on each plot bounds the region where the disc is thin, i.e. H < r (above the curve); 
(3) long-dashed line bounds the region where where thermal gas pressure is small compared to magnetic pressure and 
turbulent stress (below the curve); (4) short-dashed line separates the regions where the disc Thomson optical depth 
t c < 1 (above the line) and r c > 1 (below the line); (5) long-dashed and dotted line separates the region where the 
disc free-free optical depth f s > 1 (above the line) and fg < 1 (below the line); (6) short-dashed and dotted line 
separates the regions where the disc effective optical depth ^TgT c > 1 (above the line) and y/fsr c < 1 (below the line); 
(7) dotted line separates the regions where the disc free-free optical depth exceeds Thomson optical depth, fs > t c , 
(above the line) from where ffr < t c , (below the line). 
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Fig. 2. Plots of conditions when our model is valid. Plots are done in the plane of Bio an d r/r g for M = 10 Mq and 
S = 5/4. All notations are the same as in Fig. ^ 



The energy transfer due to repea ted scatterings (Comptonisation process) is characterized by Compton y(u) pa- 
rameter l|Rvbicki fc Lightmari Il979|l 

V = 2Max(r cs ,r c 2 s ), (40) 

where the optical depth for Thomson scatte ring t bs (V) must be measured from the effective absorption optical thickness 
t*(v) ~ 1 and is given by (formula [7.42] in lRvbicki fc Lightmanl jl979)) 

r,(^( ^ I'' 2 - (41) 

If yiy) <C 1, photons do not change their initial frequency v in the process of repeated scatterings before they escape 
the surface of the disc. While if y > 1 the Comptonisation effects become significant. We calculated values of y 
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using expressions l|4UI) . (|41|l and (|39|) . One can see that y is always monotonically rising with frequency, therefore, 
the Comptonisation effect at higher frequencies is always more significant than at lower frequencies. On the other 
hand, there is very little radiation at hv > 5kT due to exponential cut off in thermal spectra. It turns out that in 
most cases y(v) <C 1 for optically thick disc models and for hv < 5kT. Exceptions are the cases of high accretion 
rates l^/e <; 1. In those cases, inner parts of the accretion disc (r < 30r g ) can have y > 1 and Comptonisation 
will influence the highest frequencies of the disc spectrum. Ignoring these exceptions, we did not take into account 
Comptonisation in the following calculations and assume coherent scattering. This assumption has been checked a 
posteriori for self-consistency. 

In the case of coherent scattering the approximate expression of the radiative flux per unit surface of an optically 
thick medium is given bv llRvbicki fc Lightmanl ll979l) 

p v ~ I (42) 

( e W*r _ i) (i + + Kt / Ks („, T)) ' 

F v approaches black body spectrum -kB v (T) in the limit Kg ^> kt and becomes modified black body spectrum 

F v = 2itB 1/ (T)^k s / Kt (43) 

in the limit Kg <C kt- A part of the disc can emit black body spectrum at lower frequencies v <C vq and modified 
black body spectrum at higher frequencies v 3> vq, where vq = vq(p,T) is defined such that ks(vq) = kt- However, 
if frequency averaged free- free opacity is larger than Thomson opacity, Rg{T,p) > kt, then hvo > kT and almost all 
the radiation is emitted as a black body spectrum. In the opposite limit, Kg(T,p) <C kt, one has hvo -C kT and the 
spectrum is mostly modified black body, transiting to black body only at lowest frequencies, v < vq. Dotted line on 
Figs. QH2I separates regions with Kg > kt and Kg < kt- We see that the optically thick disc has Kg <C kt and emits 
a modified black body spectrum in the inner parts but may become absorption dominated in the cooler outer parts, 
where black body spectrum will be emitted. At lower accretion rates, Ie/(2€) < 1CP 4 , all the surface of the disc will 
emit black body spectrum (with different T at different radii, of course). The locations of the regions with mostly 
black body and modified black body spectra over the disc radii are fairly insensitive to the black hole mass M. The 
surface temperature T in formula (|42|l should be determined by equating the total emitted flux, the integral of F v 
over all frequencies, to the half of the heat production rate Q per unit disc surface (half is to account for two surfaces 

of the disc). Introducing variable x — —- one can write this energy balance condition as 



kT 



Q 4**£r ^ . (44) 

'° (e x - 1) (l + y/l + K T /Kg ) 



2 h 3 c 2 



2TT 5 k 4 

Further introducing Stefan-Boltzmann constant a = ^— -r into right hand side of equation (14411 and effective 

temperature Q/2 = aT^ s (equation l|26ll ) into left hand side of equation (|44|l we can transform equation (|44|l to 

T 4 tt 4 /15 



(45) 



o 1 + \j\ + kt/kh • 



7T 4 



Since — = / , T is always larger than T c g. This is in accordance with general thermodynamic argument that 

15 Jo eX - 1 

the black body is the most efficient emitter of all. Equation l|45f) together with expression l|39fl for Kg, expression l|2U|l 
for p, expression l|26|) for T e g, and kt = 0.4cm 2 g~ 1 has been solved numerically to determine T(r). After one knows 
T = T(r), it is possible to integrate F v (r) l|4^|) over the disc surface to obtain the spectral distribution of the total 
energy emitted by the disc 

E v = 2 [ 2irrF u dr, (46) 

^ ^in 

where the 2 accounts for the two surfaces of the disc. 

When a significant interval of radii exists where the emitted spectrum is a modified black body, e.g. hvo < kT, it 
is possible to get an approximate analytic expression for E v . For v > vq we use expression ()43|l for F v , which becomes 

F v = 2.6 x 1CT 3 CTg T5/y/^/ 2 ■ e ~ X (47) 
scm 2 Hz ■y/l - e ~ x 
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Integrating expression (|47|l over frequencies by integrating from < x < oo, we obtain 

^ =8 .2x 10 7 ^T 9 /V /2 , 
2 s cm 2 

which can be solved for the temperature using expression l|12|l for Q. Let us denote the temperature found in this way 
by T s . The expression for T s is 

// \ 8/9 / R \ ~ 4 / 3 

= ^"(wb) * 

-8/3+4(5/3 

g 8/9 . (48) 



10r„ 



Now we use T s to substitute in the expression l|47|) for F v , and then to expression l|46|) to obtain E u . It is convenient to 

hv dT 

switch from the integration in r to the integration in x in equation l|46|l . We do so by writing dx = — - — - — — dr and 

fcT~ dr 

expressing x through r for a given i/ from equation l|48|l . This procedure can be done analytically if one puts Q = 1, 
that is our analytical approximation do not describe spectrum emerging close to the inner edge of the disc, where Q 
deviates from 1 significantly. Carrying out calculations one obtains 

4-35 

n.73 erg /Z E \ 3p ^ 
E v = 4.2 x 10 44 x 10" — f -F x 

sHz \2e J 

( Bm\ ^ x (49 ) 

8 \10 4 GJ (2-S) v ' 

3(25-3) £~ x 

x 4 < a - 2 > , dx, 



where x- m = - — t — - and x ou t = , — , -. We take for typical estimates r m = 10r g and r out = 1000r g as the 



hv hv 

kT s (r in) ° Ut kT s (r out ) 

inner and outer edges of the disc. Although inner parts of the disc contribute significantly to the total emitted power 
and determine the most energetic part of the spectrum, the calculation of the spectrum there must be performed by 
taking into account factor Q, not to mention relativistic effects. The outer extension of the disc at 1000r g is somewhat 
arbitrary, but the disc beyond 1000r g is too cool to be described by our simple radiation model and, in the case of 
AGNs, even the existence of the disc for r > 1000r g is questionable because of the instability to the gravitational 
fragmentation. The value of the integral in formula 149fl decreases exponentially for x ln > 1. This corresponds to an 
exponential cutoff in the spectrum for hv > kT(r- m ). If x- m <ti 1 but x out > 1, then the value of the integral is almost 
independent on v and is a slowly varying function of 5. We see that E v cx ^( 4<5 ~ 5 )/( 8 ~ 415 ) in this case. Thus, E v is rising 
for 5 > 5/4 and declining for 6 < 5/4. If both x- m <C 1 and x out <C 1, then it is possible to see from the expansion of the 
integral in expression (|49|l that E v cx v. At frequencies below hv = kT(r out ) the whole disc surface would contribute 
with the low frequency tails of modified black body spectra, which scale as F v oc v (equation (|47|) ~). Therefore, it is easy 
to understand the scaling E v cx v for hv < kT(r out ). However, we do not see the latter spectral index in calculated 
spectra because Kg becomes comparable to kt already at the frequencies where x out > 1. The E v cx I /( 4<5_5 )/( 8 ^ 4 ' 5 ) 

law extends down to the frequency at which a;; n = XQ S {r- ln ), where xq s = ^— is found by equating Kg — kt- For 

kT s (r) 

xo s <C 1 one obtains using expression 139|) for Kg together with expression (|48|l for T s and expression l|20|l for p 

/ / \ -46/27 , / n \ 32 / 9 

Zos^l.SxltnM^) Af 46/27 ' Bl ° 



Vl0r g 



2eJ 8 V 1 4 G / 

46/9-325/9 

g- 46 / 27 . (50) 



For x < Xq s spectrum gradually transits to the sum of a local bla ck body with T = T e g, which has E v cx v 1 ^ 3 
( Sha kura fc Sunvaevl Il973 IShapiro fc Tcukolskv l Il98a iKrolikL Il999|) . Finally, for hv <C kT(r out ) the spectrum is the 
sum of cx v 2 low frequency black bodies of v < vq. When the outer part of the disc, where Kg > kt, is sufficiently 
truncated then the E v cx z^ 1 / 3 part of the spectrum may be absent and the spectrum will transit directly from 
E v cx j,( 4<5 - 5 )/( 8 - 415 ) power to E v cx v 2 power. 

In summary, we see that magnetically dominated accretion discs have power law spectra with the spectral index 
depending on the radial distribution of magnetic field strength such that, E v cx ;y( 4 ' 5 ~ 5 )/( 8 ~ 4 ' 5 ). This contrasts the 
standard weakly magnetized a-disc which shows a declining modified black body formed from the inner radiation 
dominated disc with E v cx v~ 2 / 5 . 
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4.2. Modified black body spectrum in a standard disc 

vd In E 

As a side remark we note that the value for the spectral index 7 = — close to found for the latter regime of 

accretion disc bv lShakura fc Sunvaevl l)l973|) (text on pag e 349 after equation [ 3.11] o f that work) is different from ours 
7 = —2/5. It is easy to follow the exact prescription of IShakura fc Sunvaevl l)l973|) . namely calculate integral [3.10] 
in their work for spect r um [3 .2] and temperature [3.7]. As a result we obtain 7 = —2/5 rather than 7 = 0.04 given 
in lShakura fc Sunvaevl lll973l). We need to point out this dis crepancy because it is widely stated in many text books 
on accretion discs ijShapiro fc Teukolskvl Il983t iKrolikl Il999|) with the reference to IShakura fc Sunvaevl l|l973l) that 
high luminosity accretion discs has almost flat plateau in its spectrum before the exponential cut off corresponding to 
kT(r s ). However, the flat spectrum E v cx v 1 / 29 is produced by part (b) of the standard disc mode l, where gas pressure 
domin ates over radiation pressure. The spectral index 7 = 1/29 is close to the 7 = 0.04 given in lShakura fc Sunvaevl 
lll973h but the radial dependence of the surface temperature in zone (b) is T s cx r 29 / 30 rather than T s cx r 5 / 3 given by 
their formula [3.7]. Thus, the standard a-disc possessing both (b) and (a) zones should have spectrum steepening from 
plateau to cx u~ 2 / b and then exponentially cutting off at the temperature of the inner edge. Because the intervals of r, 
where approximate analytic expressions for emitted spectrum are valid, do not typically exceed one order of magnitude 
(the same is true for our disc model as well), one does not see "pure" extended power laws when c alculating spectra 
numerically by using general expressions l|42|l and l|44l) . For example. IWandel fc Petrosianl {l988) found 7 = —2/5 
slope in the narrow interval between 1000 A and 1450 A by numerically integrating disc spectra. 



4.3. Results of spectrum calculations 

We present results of the simplified analytical integration of the spectrum using equation (|47|l as well as more exact 
numerical integration using equation l|42|l . solving for T from equation l|45|) and integrating equation (|46|l . Function 
Q{r) was kept in numerical calculations, so the results are applicable to the innermost parts of the disc, where the 
most of energy is radiated. For a given M 8 and Ie/c, an optically thick magnetically dominated discs exist within 
5r g < r < 1000r g only for S in the interval of about 1 to 1.4. In Figs. ISTfTTlwe illustrate models for the four choices of 
parameter sets: 

1. M s = 1, f| = 0.1, 8 = 5/4, B w = 3 x 10 3 G, 

2. M 8 = 1, l f e = 0.1, 5 = 1, Bio = 3 X 10 3 G, 

3. M 8 = 1, §| = 0.1, 6 = 1.4, B w = 5 x 10 3 G, 

4. M s = 1, t| = 10~ 3 , 8 = 5/4, B w = 7 x 10 2 G. 

The dependencies of H on r given by equation i|18fl for four parameter sets are plotted in Fig. [31 For comparison we 
also plot the half-thickness H in the Shakura-Sunyaev model of the disc with the same accretion rate (parametrized by 
Ie/c) aud the same mass of the central object. We use approximate analytic expressions for the parameters of the disc 
(H , p, E, T mnr {) in the radiation dom inated zone (a) of the Shakura-Sunyaev model and thermal pressure dominated 
zone (b) llShakura fc Sunvaevl Il97l . The magnetically dominated disc is thicker than the standard disc. For higher 
accretion rates, the standard disc has a concave shape due to the transition from inner zone (a) to intermediate 
zone (b), which allows illumination of the outer parts of the disc by the inner hotter and brighter parts. Magnetically 
dominated disc has convex shape, which exclude such illumination. 

In Fig. 2] we plot the dependencies of the column thickness through the disc E (equation j!9Jl) on the radius and 
also compare to Shakura-Sunyaev standard disc. The magnetically dominated disc is much less massive than the 
standard disc. Both E and p are smaller for magnetically dominated discs, and only in the inner ~ 10r g are the 
densities comparable. 

The dependencies of mid-plane temperature T mp( j on radius given by equation Ij28(l are shown in Fig. 03 On the 
same figure we also plot T mpc i in the Shakura-Sunyaev model and T c g given by equation i|26|) . which is the same for 
magnetically dominated and standard discs. Because of the lower column density of the magnetically dominated disc, 
T mp d is less than for standard cv-discs. 

The dependencies of magnetic plus turbulent pressure, ~ B 2 /(47r), radiation pressure in the disc mid-plane P r ad = 
aX^pj/3, and thermal pressure in the disc mid-plane Pth = ufcT mp( i are presented in Fig.0 We see that the assumption 
of magnetic pressure dominance is well satisfied for our models except in the innermost regions, r ^ 10r g , for higher 
accretion rates Ze/(2c) = 0.1, where radiation pressure becomes comparable to the magnetic pressure. The latter fact 
limits the existence of magnetically dominated regime in the innermost parts of accretion discs for higher luminosities. 
Plasma parameter (3 defined as f3 — 87r(P rac i + P t h)/-B 2 decreases with radius and varies from ~ 1 to ~ 10~ 2 in our 
models. 

In the limit (3=1 magnetic pressure is comparable to the largest of radiation or thermal pressures and our strongly 
magnetized disc model transforms into Shakura-Sunyaev model with a = 1. If 6 = 3/4 in our model, then the radial 



V.I. Pariev, et al.: A strongly magnetized accretion disc model 



15 



\ 

o 







-2 

2 







III 


I I I I 


i i i i i I 




~ B 10 = 3xl0 3 G 








l E /(2e) = 0.1 








/ 5 = 5/4 






1 


- 

/ 

/ 




/ / — 

y / 





— //- -~~ 






-1 


- ' 






-2 






I | I I I | ^ 


o 
c. 




I I 




1 1 .5 


2 


2.5 3: 




B 10 = 3xl0 3 G 








1 l E /(2e) = 0.1 






1 


- 5 = 1 






- 


// 

UTl / 


— 











-1 


- 

r i i i 1 i i i i lii 


ill 


i i i i i I 


-2 


1 1.5 


2 


2.5 3 




log 


(r/rj 







i i i i i i i i i r 

B 10 = 7xl0 2 G 
l E /(2e) = lE-3 
5 = 5/4 




i i i 
1 



1.5 



i i i 
2 



2.5 



log(r/r 



Fig. 3. Dependencies of the half-thickness of the disc H on radius for magnetically dominated disc (solid line) and 
Shakura-Sunyaev disc with the same Ie/ 'e and Ms parameters and viscosity parameter a = 0.1 (short-dashed line) and 
a = 0.01 (long-dashed line). The breaks in the curves for the Shakura-Sunyaev disc occur at the interface of zones (a) 
and (b) and are the results of using approximate analytic expressions in zone (a) and zone (b). For l-^/{2e) = 10~ 3 
zone (b) extends down to the inner edge of the disc. 



scaling of the magnetic and turbulent pressures, -B 2 /(47r), is the same as that of the radiation pressure inside the 
disc, etX^pj, aside from the factor Q. If S = 51/40, then the radial scaling of B 2 /(4tt) is the same as that of the 
thermal pressure inside the disc, nfcT mp d. Therefore, by choosing (5 = 3/4 and adjusting the magnitude of Bio one can 
construct the model with approximately constant f3 in the zone where radiation pressure exceeds thermal pressure. 
By choosing S = 51/40 one can construct constant /3 model in the zone where thermal pressure dominates radiation 
pressure. We illustrate this in Fig. [7| where we show the dependencies of pressures, p, H, and S on r for our model 
with 8 = 3/4, Ze/(2c) = 0.5, and M% = 1, and for the Shakura-Sunyaev model with a — 1 for the same accretion rate 
Ze/(2e) and M 8 . The transition from zone (a) to zone (b) in this Shakura-Sunyaev model occurs at r a b = 360r g . The 
breaks on the curves corresponding to the Shakura-Sunyaev model occur at r = r a b in Fig0 We adjusted B w such 
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Fig. 4. Dependencies of the surface density E on radius for magnetically dominated disc (solid line) and Shakura- 
Sunyaev disc with the same Ie/ 'e and parameters and viscosity parameter a = 0.1 (short-dashed line) and a = 0.01 
(long-dashed line) . The breaks in the curves for the Shakura-Sunyaev disc occur at the interface of zones (a) and (b) 
and are the results of using approximate analytic expressions in zone (a) and zone (b). For Ze/(2c) = 10~ 3 zone (b) 
extends down to the inner edge of the disc. 



that the magnetic pressure will be in equipartition with the radiation pressure in our model. Then, as it is seen from 
the top-left plot in Fig. thermal pressure is less than magnetic and radiation pressures for r less than some r c and 
exceeds magnetic and radiation pressures for r > r c , so our model is not applicable for r > r c . The subsequent three 
plots show that r c sa r a b. Two right plots and bottom-left plot in Fig.0show that p, H, and E in our model for r < r c 
are very close to p, H, and E in Shakura-Sunyaev a = 1 model in the radiation pressure dominated zone r < r a b. 
A similar conclusion holds for the transition of our model with 6 = 51/40 to a Shakura-Sunyaev zone (b) model for 
r > r a b and (3=1. The radiation spectra of our model in the limit (3=1 also approach that of the Shakura-Sunyaev 
model, as shown by direct numerical calculations. The power law modified black body spectrum E v oc j/( 4,5 ~ 5 )/( 8 ~ 415 ) 



V.I. Pariev, et al.: A strongly magnetized accretion disc model 



17 



~i — r~i — n — r~n — n — r~ n — n — r~n — n — r 




l E /(2e) = 0.1 

5 = 5/4 



I I I | I I I 
1.5 



I I I | I I I 
2.5 



3~ 




= 3xl0 3 G 
l E /(2e) = 0.1 
6=1 

i i I i i i i I i i i 



i i i i i i i 



1 1.5 



2.5 



log(r/r ) 



~th — i i i — r~r 




B 10 = 7xl0 2 G 
l E /(2e) = lE-3 
5 = 5/4 



1.5 



2.5 IT 




B 10 = 5xl0 3 G 
_l E /(2e) = 0.1 
5=1.4 



1.5 



2.5 



log(r/r ) 



Fig. 5. Dependencies of temperatures on radius: T mp d for magnetically dominated disc (solid line); T m pd for Shakura- 
Sunyaev disc with the same Ze/c and Mg parameters and viscosity parameter a = 0.1 (short-dashed line) and a = 0.01 
(long-dashed line); T e g (dashed-dotted line). The breaks in the curves for the Shakura-Sunyaev disc occur at the 
interface of zones (a) and (b) and are the results of using approximate analytic expressions in zone (a) and zone (b). 
For Ze/(2c) = 10~ 3 zone (b) extends down to the inner edge of the disc. 



derived in section FTTI becomes E v oc v~ 

-2/5 for 5 = 3/4 an d E v cx v 1 ' 29 for S = 51/40, which is coincident with the 
modified black body power laws for the Shakura-Sunyaev zone (a) and zone (b) spectra (see section l4~2t . 

The dependencies of optical depths through the half disc thickness on radius are shown in Fig.IHlfor four parameter 
sets. The three curves plotted are: fg given by equation l|35|) . t c = ktS/2, and the effective optical thickness f* = 

-z\J («t + Sff)Sff. For S — 5/4, the effective optical thickness f* is almost constant throughout the disc, but when 8 

deviates from 5/4, f* starts to approach 1 either at the inner or at the outer edge of the disc and so our model breaks 
down at those radii. With the decrease of the accretion rate, the disc becomes cooler and denser so the absorbing 
opacity rises and becomes larger than the scattering opacity in the outer parts of the disc. 
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Fig. 6. Dependencies of pressures on radius for parameters of disc considered in section Magnetic plus turbulent 
pressure B 2 /(4tt) is plotted with a solid line, P ra d is plotted with a short-dashed line, and Pth is plotted with a 
long-dashed line. 



In Fig.Elwe show: T c s(r) given by equation T mp d(r) given by equation T s (r) given by equation and 
T(r) by solving equation (|45J) numerically. Note that T e s and surface temperature T are always smaller than the T mp d 
for an optically thick disc. For low accretion rate, Ze/(2c) = 10 -3 , Rg > kt and T w T c g. In this case, T s is ill defined 
and values of T s < T c g are unphysical on the plot for Ze/(2c) = 10 -3 and also for r > 100r g on the plot for 6 = 1, 
Ze/(2c) =0.1. The temperature T s becomes a good approximation for T when T s > T e g (scattering dominates over 
absorption in the surface layer). Unlike the values and slope of T mp d(r), which substantially increases with increasing 
5, the value of T is less sensitive to 5: only inner parts of the disc becomes slightly hotter for larger values of 8. Both 
T and T e g are changed significantly when the accretion rate or mass M are changed. 

Fig. 1101 shows the results of calculating Comptonisation y parameter according to equations H4l)|) and 1)41(1. We 
conservatively set x = 5 for the calculation of y. Then, y is the function of radius alone. On the same figure we also 
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Fig. 7. Comparison of /3 = 1 limit of magnetically dominated disc model with 8 = 3/4, l^/(2e) — 0.5, and Mg = 1 
and Shakura-Sunyaev model. Top-left plot shows the radial profiles of magnetic plus turbulent pressure = B 2 /(4tt) 
(solid line), P ra( j (short-dashed line), and Pth (long-dashed line) in the magnetically dominated model. Bottom- left, 
top- right, and bottom-right plots show the radial profiles of p, S, and H in the magnetically dominated model (solid 
lines) and the Shakura-Sunyaev model with a = 1 (dashed lines). The parameters l-&/{2e) and M§ are the same for 
magnetically dominated model and Shakura-Sunyaev model. 



show y in the regime of modified black body spectrum, using T s and writing the simplified version of equation l|40l) as 
AkT kt 

y = — . We see that Comptonisation is not important for our models even in the inner disc. 

m c cr Kg 

Energy spectra E v are presented in Fig. ^2 We normalized frequency to the characteristic frequency of an effective 
black body from the inner disc, namely, we plot versus v/v e Si where 
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Fig. 8. Dependencies of optical depth through the half disc thickness on radius for parameters of disc considered in 
section Dashed line is r c , dashed-dotted line is fg, and solid line is f*. 



so the spectra plotted cover the range from EUV to infrared for M = 10 8 Mq. We checked that the total thermal 
energy emitted from the disc between r in and r out calculated as an integral of the spectrum over frequencies is equal 
to the surface integral of the dissipation Q (expression 1121) ): 



E = E v dv = 



2irrQ dr 



9.5 x 10 45 erg s" 1 — M 8 — x 

e rir, 



(51) 



V.I. Pariev, et al.: A strongly magnetized accretion disc model 



21 



l l | l l l l | l l l l | l l l l | l l l l 




log(r/r g ) log(r/r g ) 

Fig. 9. Dependencies of temperatures on radius for parameters of disc considered in section^] Short dashed line is T e g, 
long dashed line is T mp d, dotted line is T s , and solid line is T. The latter two temperatures are the surface temperatures 
of the disc calculated in section 0] 



All spectra shown in Fig. 1111 were computed by integrating from r- ln — 3.1r g to r out = 1000r g with the factor G(r) 
taken into account. We also show the spectrum calculated by using the approximate analytic expression l|49|l for a 
modified black body in its regime of validity (x; n > XQ S (ri n )). Because the analytic expression was obtained by setting 
Q = 1 in equation i|48|) for T s , it overestimates the temperature in the inner parts of the disc by a factor of a few and 
does not describe the high energy part of the spectra correctly. The lower frequency at which the sum of modified 
black bodies is still a good approximation, increases with the overall increase of absorption in the disc. For parameter 
set 01 above, with Ie/(2c) = 10~ 3 , pure modified black body cannot be used at all, so the corresponding panel in 
Fig- IHI does not show a second curve. The spectrum for 5 = 5/4 (solid line) shows flat plateau extending by more 
than an order of magnitude in frequency in accordance with analytical result. Although the declining top part of the 
spectrum for (5 = 1 and rising top part for S — 1.4 are apparent, the interval of frequencies, where a modified black 
body approximation works, becomes small and blends with the cx v 1 ^ 3 spectrum of the sum of local black bodies. 
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Fig. 10. Dependencies of Comptonisation parameter y on radius for parameters of disc considered in section^] Solid 
line is for y calculated using temperature T, exact expression <|41ll for t cs , and hv = 5kT(r). Dashed line is for y 
calculated using T s , assuming that r cs ~ (tt /ts) 1 / 2 , and hv = 5kT s . 



Thus no dependence on S is evident here. The top right plot for low luminosity Ze/(2e) = 10~ 3 is indistinguishable 
from the sum of the local black body spectra. All spectra behave like cx v 2 for low frequencies. In Fig.^2 we also show 
spectra of a Shakura-Sunyaev a-disc with the same Ze/c an d Mg parameters to compare with plots of a magnetically 
dominated disc with viscosity parameter a — 0.1. These spectra were calculated in the same way we calculated the 
spectra of the magnetically dominated disc: first we found the surface temperature T(r) by solving equation l|45JI with 
the p(r) profile taken from standard a-disc model, and then integrated equation (|46|l with F v given by expression (|42|) . 
For low accretion rates of order of 10 -2 of Eddington accretion rate and s maller, the spectra of Shaku ra-Sunyaev 
disc are very close to the sum of local black bodies with temperatures T e g(r) IjWandel fc Petrosianl Il988» . In general 
the spectra of our magnetically dominated discs are close to the spectra of Shakura-Sunyaev discs, so it seems to be 
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Fig. 11. Spectral energy distribution for the total flux from the disc. Frequency is plotted in units of u e g — 
kT c s(10r g )/h. Exact values of E u calculated using temperature T are plotted with solid lines; values of E v calcu- 
lated from analytical expression l|49(l are plotted with short-dashed lines. The latter are shown only for frequencies at 
which x m is larger than the minimal value of xq s (t), which is achieved at about 5r g to 10r g . Spectra of Shakura-Sunyaev 
discs are plotted for viscosity parameter a = 0.1 with long-dashed lines. 



difficult to distinguish between them observationally. This means that sources with discs previously thought to be 
thermally supported could actually be magnetically supported. 



5. Discussion and Conclusions 

We have found self-consistent solutions for thin, magnetically supported turbulent accretion discs assuming the tan- 
B 2 

gential stress — a(r)- — . When compared to the standard a-disc models (Sh akura k. Sunvaevl ll9T3^l magnetically 

47T 

dominated discs have lower surface and volume densities at the same accretion rate. This is due to the more efficient 
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angular momentum transport by supersonic turbulence and strong magnetic fields than the subsonic thermal turbu- 
lence of the standard model. As a result, magnetically dominated discs are lighter and are not subject to self-gravity 
instability. In the limit of plasma (3 = 1, magnetic pressure is comparable to the largest of radiation or thermal 
pressures and our strongly magnetized disc model transforms into the Shakura-Sunyaev model with a — 1 . 

When we derived the disc structure, we made no explicit distinction between turbulent and magnetic pressure 
support and angular momentum transfer. As such, our model would be valid in any situation in which the magnetic 
and turbulent kinetic energies are comparable to, or greater than the thermal energy density. The assumption that the 
kinetic and magnetic energies are nearly comparable is natural because turbulence should result in the amplification 
of small scale magnetic fields in highly conducting medium due to dynamo action. Typically, in a sheared system, 
the magnetic energy can be even slightly larger than turbulent kinetic energy since the magnetic energy gains from 
the additional shear. We find that the thermal spectrum from the surface of the magnetically dominated disc in the 
optically thick regime is close to the spectrum of the standard Shakura-Sunyaev disc. 

The issue arises as to how the magnetic field could reach sonic or supersonic energy densities. To obtain sonic 
turbulence and produce a (3 = 1 disc, the MM might be sufficient. To obtain a [3 < 1 supersonic turbulence may require 
something else. One possibility in AGN appeals to the high density of stars in the central stellar cluster surrounding 
AGN accretion discs. Passages of stars through the disc might be an external source of supersonic turbulence analogous 
to the supernovae explosions being the source of supersonic turbulence in the Galaxy. Stars pass through the disc with 
the velocities of order of Keplerian velocity, which is much larger than the sound speed in the disc. We consider the 
support of turbulence by star-disc collisions in Appendix B and find that statistically speaking, star-disc collisions 
are unlikely to provide enough energy to sustain supersonic turbulence in most AGN accretion discs, however the 
possibility remains that a small number out of a large population could become magnetically dominated. 

Indeed whether a disc could ever really attain a magnetically dominated state is important to understand. The 
present answer from simulations is not encouraging, but not completely ruled out. Further global MHD simulations 
of turbulence in vertically stratified accretion discs with realistic physical boundary conditions are needed along with 
more interpretation and analysis. Magnetic helicity conservation for example, has not been fully analyzed in global 
accretion disc simulations to date, and yet the large scale magnetic helicity can act as a sink for magnetic energy since 
magnetic helicity inverse cascades. 

As an intermediate step in assessing the viability of low (3 discs, it may be interesting to assess whether they are 
stable. One can take, as an initial condition, the stationary model of the magnetically dominated accretion disc given 
by expressions (|15fl , (|18H20|) with the initial magnetic field satisfying all constraints of our model and falling into the 
shadowed regions on plots in Figs. 

One point of note is that mag netically dominated discs may be helpful (th ough perhaps not essential, if large scale 
magnetic fields can be produced l)BlackmarA 120021: lBla,ckman Rr. Pa.rievi [2003)) in explaining AGN sources in which 
40% of the bolometric luminosity comes from hot X-ray coronae. If the non-thermal c ompo nent i n galactic bla c k hole 
sources is attributed to the magnetized corona above the disc (e.g.. lDi Matteo et alJ l|l999|) . also iBeloborodovl l|l999l) 
discusses possible altern atives), then magnetically dominated discs can naturally explain large fractions, up to 80 % 
ijDi Matteo et al.lll999|) . of the accretio n power being transported into coronae by magnetic field buoyancy (although 
(3 ^ 1 disc solutions are also possible, i Merlonil 12003)). Though coronae can form in systems with high j3 interiors, 
the percentage of the dissipation that goes on in the interior vs. the coronae could be (3 dependent. 

The main purpose of our study was simply to explore the consequences of making a magnetically dominated 
analogy to Shakura and Sunyaev, and filling in the parmeter regime which they did not consider. In the same way 
that we cannot provide proof that a disc can be magnetically dominated, they did not present proof that a disc must 
be turbulent, but investigated the consequences of their assumption. We also realize that the naive a disc formalism 
itself can be questioned and its ultimate validity in capturing the real physics is limited. Nevertheless it still has an 
appeal of simplicity. 

Finally, we emphasize that our model does not describe dissipation in the corona and interaction of the corona 
with the disc. Further work would be necessary to address relativistic particle acceleration and emission, illumination 
of the disc surface by X-rays produced in the corona and subsequent heating of top layers of the disc, and emergence 
of magnetized outflows. 
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Appendix A: On cyclotron emission 

Since the characteristic temperature inside the disc (T mp( j given by equation (|28[l ) is non-relativistic, cyclotron emission 
of an electron occurs at frequencies close to multiples of the gyrofrequency u> = swg, where s = 1,2,3,..., and 

However, because the magnitude of the magnetic field varies across the disc, the resulting emission blends many 
discrete gyrolines. The typical range of cyclotron emission is cm radio waves for AGNs and submillimetre to infrared 
for stellar mass black holes. Characteristic plasma parameters in our disc for the case of supermassive black holes 
are similar to those encoun tered in solar chrom osphere, where plasma effects are important for the generation and 
propagation of radio waves IjZheleznvakovi 119701) . One should expect that collective plasma effects will influence the 
cyclotron radiative process at such low frequencies. The plasma frequency is 
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and the ratio of cyclotron to plasma frequencies is 

/ \ -3+25 

We see that typically lob ~ wl 5 so that cyclotron emitted radiation can not propa gate for some disc p arameters. 
However, even if lob 3> wl, the plasma affects cyclotron radiation. As summarized in lZhelezniakovl (^996) collective 
effects suppress the emission on the first harmonic, s — 1, such that it becomes of order of the emission on the second 
harmonic, s — 2. The emissivity on higher harmonics, s > 2, is smaller by the factor (kT/m e c 2 ) s ~ 2 than on the 

B 2 3/cT 

47T TO;? C 4 

c 2 VAC s m p /m c , where c s ~ 3fcT/m p . In vacuum, note that the emissivity on the first harmonic is (m e c 2 /kT) times 
larger than the emissivity on second harmonic. The latter condition is narrowly satisfied for small radii of our optically 
thick disc models and the larger the r the better it is satisfied. 

Cyclotron self-absorption also occurs in narrow lines centred on multiples of u>b- At some frequency w, emission and 
absorption occurs only in spatially narrow resonant layers inside the disc, where the magnetic field strength matches 
the frequency, i.e. slob(B) — uj is small. The width of emission and absorption frequency intervals is determined mainly 
by the thermal Doppler shifts Auj/w ~ \/kT/m P c 2 l(Zhelezniakovi fl99fil) . The width of such resonant layers can be 
estimated as ~ Hy/ kT/m c c 2 . IZhelezniakovl l|l99fil) (chapter 6) gives the expression for the optical thickness through 
such gyroresonance layers on the second harmonic uj — 2lub, which for our disc is 



second and first harmonics. This occurs only at high enough plasma densities n 3> ^— r, which translates into 
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For all our models T cyc is always very large and the emission is always strongly self-absorbed. Cyclotron photons are 
not subject to Compton scattering by free electrons, since the wavelength of the emission is always larger than Debeye 
radius in plasma, so electrostatic shielding of charges prevents them from scattering. Under such circumstances the 
resulting cyclotron flux from each gyroresonance layer is that of the black body with the local plasma temperature 
in the gyroresonance layer. Due to the overlapping of all layers, the resulting spectrum is a black body spectrum of 
width ~ 2u>b- Since Tiuob <C kT, the total flux of cyclotron emission from the disc surface is negligibly small. 



Appendix B: Star-disc collisions as possible source of turbulence 



When a star passes through a disc, it creates strong cylindrical shock propagating in the surrounding gas in the disc. 
The aftershock gas is heated to temperatures exceeding the equilibrium temperature in the accretion disc. As the 
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shock weakens, this heating decreases until at some distance from the impact point the incremental heating becomes 
comparable to the equilibrium heat content. The scale substantially affected by a star passage is x ~ R*v^/ 'c s , much 
larger than the radius of the star . The shock front can become unstable and turbulence can occur in the aftershock 
gas. The heated gas becomes buoyant, rises above the disc and falls back because of gravity. Fall-back occurs with 
supersonic velocities and can further excite turbulence. Turbulence will derive energy from both heating by star 
passages and shear of the flow. The energy, which can be derived from shear, is equal to Q given by expression (|12|) . 
It is possible that star-disc collisions might mainly be a trigger for the available shear energy to be converted into 
supersonic turbulence, and additional energy deposited into the disc by star-disc collisions is negligible. However, it 
seems unreasonable that the star-disc collisions can influence the global structure of the accretion disc unless the 
energy deposition from them is some fraction of the energy necessary to sustain turbulence level Q in the disc. 
The energy deposition rate by stars per unit surface of the disc is 

Q* « nRl'Ev^n^v^, (B-l) 

where w vk is the typical velocity of stars at radial distance r, R* sa Rq is the average radius of stars, n* = n*(r) is 
the number density of stars, £ is the surface density of the disc given by expression i|19|) in our model. Only accretion 
discs orbiting supermassive black holes M ^ 10 6 Mq can be influenced by star-disc collisions. 

The resolution of observations is only enough to estimate the number density of stars at about 1 pc for M32 
and M31 and about 10 pc for nearest ellipticals. In line with these o bservations we ass ume a star density n(l pc) « 
10 4 — 10 6 Mqpc -3 at 1 pc distance from the central massive black hole ijLauer et al.lll995|) . To estimate for r < 10 3 r g 
we need to rely on the theory of cent ral star cluster evolutio n. The gravitational potential inside the central 1 pc will be 
always dominated by the black hole. lBahcall fc Wold 1)1976|) showed that, if the evolution of a star cluster is dominated 
by relaxation, the effect of a central Newtonian point mass on an isotropic cluster would be to create a density profile 
n oc r -7 / 4 . However, for small radii (« 0.1-1 pc) the physical collisions of stars dominate two-body relaxations. Also, 
regions near the black hole will be devoid of stars due to tidal disruption and capture by the black hole. Numerical 
simulations of the evolution of the central star cluster, taking into account star- star col l isions , star-star gra vitational 
interactions, tidal disruptions and relativistic effects were recently performed by iRauchl ((1999). Raild] l)l999|) showed 
that star-star collisions lead to the formation of a plateau in st ellar density fo r small r because of the large rates of 
destruction by collisions. We adopt the results of model 4 from IRauchl l(l999l) as our fiducial model. This model was 
calculated for all stars having initially one solar mass. The collisional evolution is close to a stationary state, when 
the combined losses of stars due to collisions, ejection, tidal disruptions and capture by the black hole are balanced 
by the replenishment of stars as a result of two-body relaxation in the outer region with a n cx r -7 / 4 density profile. 
Taking into account the order of magnitude uncertainty in the the observed star density at 1 pc, the fact that model 4 
has not quite reached a stationary state can be accepted for order of magnitude estimates. For M = 10 8 M§ Mq we 
approximate the density profile of model 4 as 

n = n 5 x 10 5 — ^ for r > 10" 2 pc, 

pc J \ lpc / 

n = ri 5 x 3 x 10 8 — ^ for 5r t < r < 10~ 2 pc, (B.2) 
pc J 

n decreases for r < 5r t , 

where r t = 2.1 x 10~ 4 pc x M^ 3 — 20r g Mg~ 2 ^ 3 is the tidal disruption radius for a solar mass star, and 715 = 
nil pc) 

n . — ; 5-. The region r < rt is completely devoid of stars. 

10° Mq/pc _j 

We see that star-disc collisions cannot excite turbulence and strong magnetic fields in the very inner part of the 
accretion disc, for r < r t , and such excitation should be weak for r t < r < 5r t . The relative width of the star depleted 
region, 5r t /r g , decreases with increasing M. For M«2x 10 9 Mq r t = 3r g and star-disc collisions happen all over the 
disc. For M < 3 x 10 6 Mq 5r t > 10 3 r g and for M < 3 x 10 5 Mq r t > 10 3 r g and star-disc collisions are unimportant 
for the structure of the accretion disc. Let us assume that it should be = fQ where the fraction / is less than 
unity but not much less than unity. Further, we use expression (|12f> for Q, expression 119|) for S, and the value for 
n in the constant density core of the star cluster, second raw in expression l|B.2f) . to substitute into Q* — fQ- Since 
the relation Q* = fQ should be satisfied for all values of r, the value of S is determined and turns out to be 6 = 3/2. 
Solving the rest of the equation for the magnitude of magnetic field Bio at 5 — 3/2 we obtain 



S i2_ = r k v in* fh ] w- 3 /^ -1/4 f 1/4 
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This value of B required by energy input from star-disc collisions should fall into the range of constraints for the -Bio 
listed in the end of section [21 

We explored all feasible range of parameters M§, Isfe, f > 10 -3 , n$ < 10 3 and found that the magnetic field 
calculated from expression l|B.3(l is always too strong to fall in the allowed range of parameters discussed at the end of 
section[3] In particular, the constraint that magnetic and turbulent pressure dominate thermal and radiation pressure 
is violated. The minimum number density of stars necessary to satisfy this constraint at the most favourable values 
of other parameters still plausible for some AGN (M§ = 40, Ze/c = 10~ 10 , / = 10~ 3 ), turns out to be « 10 6 . Such 
a high number density of stars would imply total mass in stars of order of 10 n MQ inside the central parsec from the 
central black hole. This mass exceeds observational and theoretical limitations. 



